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Abstract 

Wc discuss spherically symmetric black holes in the modified self-dual theory of grav- 
ity recently studied by Krasnov. obtained adding a Weyl-curvature dependent 'cosmological 
term' to the Plebanski lagrangian for general relativity. This type of modified gravity admits 
two different types of singularities: one is a true singularity for the theory where the fun- 
damental fields of the theory, as well as the (auxiliary) spacetime metric, become singular, 
and the other one is a milder "non-metric singularity" where the metric description of the 
spacetime breaks down but the fundamental fields themselves are regular. We first gener- 
alise this modified self-dual gravity to include Maxwell's field and then study basic features 
of spherically symmetric, charged black holes, with particular focus on whether these two 
types of singularities are hidden or naked. We restrict our attention to minimal forms of 
the modification, and find that the theory exhibits 'screening' effects of the electric charge 
(or 'anti-screening', depending upon the sign of the modification term), in the sense that it 
leads to the possibility of charging the black hole more (or less) than it would be possible 
in general relativity without exposing a naked singularity. We also find that for any (even 
arbitrarily large) value of charge, true singularities of the theory appear to be either achronal 
(non-timelike) covered by the hypersurface of a harmless non-metric singularity, or simply 
hidden inside at least one Killing horizon. 
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1 Introduction 

Ever since general relativity was formulated, various different types of modifications of the theory 
have been proposed, aiming at, e.g., obtaining better control of gravity in quantum regime or 
attempting to account for astrophysical/cosmological observations without invoking unknown, 
exotic matter fields in the cosmic inventory. Recently, Krasnov has introduced yet another 
modification of general relativity [1, 2, 3, 4, 5, 6] where the standard cosmological constant 
term is made into a function ^ of the Weyl curvature, in such a fashion as not to add extra 
degrees of freedom. The absence of extra degrees of freedom is a remarkable property which 
distinguishes the case at hand from usual modifications of gravity. The key to this result is 
to give up the metric as the fundamental field of the theory (hence the definition "non-metric 
theory of gravity" used for instance in [1, 2, 3]), and to use the formulation of general relativity 
in terms of a self-dual 2-form as the fundamental field. This idea was introduced long ago by 
Plebanski [8], and later developed by Capovilla, Dell, Jacobson and Mason [9, 10]. A metric 
can be derived from this 2-form through the imposition of suitable constraints. Krasnov's main 
result is that these constraints can be generalised to a class of theories labelled by the choice of 
all of which give dynamics for a metric with only two degrees of freedom. Such generalisation 
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had already been considered in the Hamiltonian framework by Bengtsson and Peldan [11, 12]. 
See also [13, 14, 15]. 

It is tempting to speculate that such a modification of general relativity gives an effective 
description of quantum gravity models based on the Plebanski lagrangian, e.g. the spin foam 
formalism for loop quantum gravity [18]. In this paper, however, we consider the theory only from 
a classical viewpoint, with the goal of investigating the physical consequences of the modification, 
and its viability as a theory of classical gravity. 

The nature of the modification makes gravity behave differently in regions of different Weyl 
curvature. As the Weyl curvature vanishes on conformally flat metrics, such as Minkowski or 
FLRW cosmological spacetimes, these are also exact solutions in the modified theory. The effect 
of the modification can be studied perturbing around these solutions, but one needs to go to 
second order in perturbation theory. It is then easier to investigate the classical effects of the 
modification for spherical symmetric solutions with non-zero Weyl curvature, where departures 
from general relativity are expected right away. This case has been considered by Krasnov and 
Shtanov [3]. They looked at specific choices of ^ without matter, and found solutions analogue 
of the Schwarzschild black hole.^ A characteristic feature discovered in [3] is the existence of a 
new class of singularities, where the derived metric is singular but the fundamental fields of the 
theory remain finite. 

In this paper we extend the analysis of the theory to the coupling to electromagnetism. In 
spite of the derived nature of the metric field, a natural coupling in terms of the fundamental 
2-form was found by Capovilla, Dell, Jacobson and Mason [10], and as anticipated by Krasnov 
[2], such coupling can be extended straightforwardly to the modified theory. This allows us to 
study spherically symmetric electro- vacuum black hole solutions, for which we prove staticity 
for general form of ^. We then restrict our attention to the two simplest profiles of and give 
explicit solutions. We describe how the asymptotic structure of the solutions depends on and 
for one of the two chosen profiles we find the standard Reissner-Nordst0m-de Sitter asymptotic 
structure. On such a solution, we investigate its geometry near the horizon and the way this is 
affected by the modification. 

We find that the non-metric singularities are still present and their location depends on 
the value of the charge. The number and locations of the horizons are also affected by the 
modification. With respect to general relativity, we find that more/less charge can be poured 
into the hole before the event horizon disappears, depending upon the sign of the modifying 
term: more charge when the sign of the modification is coherent with a negative cosmological 
constant, and less charge when it is coherent with a positive cosmological constant. In a sense, 
one could say that the electric charge is screened in one case, and anti-screened in the other case. 
In particular, in the first case, a sector of the theory exists where arbitrarily charged solutions 
admit a (Killing) horizon which is located inside a non-metric singularity (for this reason, it 
cannot be identified with the event horizon) and covers the central singularity at r = 0. This 
striking departure from general relativity shows that useful bounds on the physically allowed 
types of ^ can be obtained by looking at classical solutions coupled with matter. 

'^A central motivation of Krasnov and Shtanov [3, 6] is to investigate whether the modification to gravity can 
be suitable to explaining the anomalous rotational curve of spiral galaxies without appealing to dark matter, and 
they have considered profiles of $ tailored to that end. 
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The plan of the paper is as follows. In order to keep the paper as self-contained as possible, 
and also to establish our notation and conventions, the next two Sections are devoted to reviews 
of self-dual gravity theories; Plebahski's self-dual formulation of general relativity in the next 
Section and Krasnov's modified self-dual gravity in Sect. 3. We also refer the reader to the 
literature provided in the bibliography for more details on both the standard Plebahski formu- 
lation of gravity and the modification proposed by Krasnov. Section 4 describes the coupling to 
electromagnetism. In Section 5 we restrict attention to the spherical symmetric case, and show 
how to obtain the field equations. We roughly follow [3], although the notation and presentation 
differ here and there. We also reprove Birkhoff's theorem in a slightly more conventional way. In 
Section 6, we discuss the rationale for ^ and possible links to the existence of asymptotically flat 
solutions. In Section 7 we discuss the field equations for the minimal modification which allows 
asymptotic flatness, and the existence of Reissner-Nordstr0m and Nariai-Bertotti-Robinson so- 
lutions. We then focus on this minimal modification; in Section 8 we review the analytic solution 
for the vacuum black hole found by Krasnov and Shtanov in [3] , and in Section 9 we present the 
solutions for the charged case. In the final Section 10 we give an overview of the properties of 
the solutions. Throughout the paper, we set c = 47reo = 1. 



2 Self-duality and the Plebanski Action 

The notion of self-duality plays a key role in the Plebahski construction of general relativity. 
Consider a 2-form F = ^ F^i, dx^ A dx^ on a pseudo-Riemannian manifold with metric g^^ and 
non-zero determinant < 0. Introduce the identity 1 = 5pa = ^^p^^^ (here and in the following 
a square bracket on the indices means normalised anti-symmetrisation) and the Hodge star 
■k = e^'^ pa (here e is the completely antisymmetric tensor density with g'^^^s _ -|^^^ The 

latter satisfies = — 1, as can be seen immediately using 

e''%ue''\r = ^gSZ- (1) 

Then, the complex orthogonal projectors = ^(1 =b i-k) (i.e. V^"^ = , jy^jy^ = Q) define 
respectively the self-dual and antiself-dual parts of the 2-form F, 

Flip = {-P^ F),. = \ (^.^ ± ^ ^"^^ Fp^) ■ (2) 

These are respectively the positive and negative eigenvectors of i*, 

i*FW = ±F(±), (3) 

and we have F = -|- F^~\'^ An important property which will play a role in the modified 
theory is the invariance of -k under conformal transformations of the metric g^i, . 

In the following, we will work with forms valued in the Lorentz algebra so(3, 1). A typical 
example is the tetrad e^(x), defined by g^^^, = ejj^elrjjj. The tetrad is a 1-form with indices 
/ = ... 3 in the fundamental representation of the Lorentz algebra. Another well-known 

^We put the label ± of self and antiself-dual form within a bracket, to avoid confusion with SU(2) indices later 
on. In the literature one also finds the convention = ^(1 =F**), so that self and antiself are the ±i eigenvectors 
of *. 
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example is the spin connection u!j/{e), a 1-form with values in the adjoint representation of 
so(3, 1). For an object in the adjoint representation, we can define the algebraic self-dual and 
antiself-dual components, as we did in (2). This time we have to use the Hodge operator on the 
Lorentz bundle, which reads \e^'^ kl] and thus = ^((^j^x ^ ^^^"^ kl)- These two projectors 
realize explicitly the familiar homomorphism so(3, 1; C) = sur (2)©sml(2) of the Lorentz algebra 
into two su{2) algebras, which rather than self-dual and antiself-dual are more commonly dubbed 
right-handed and left-handed. To make the mapping more explicit, it is convenient to pick out 
the time direction / = 0, and define a;^^)* = w^^^'^*, with z = 1,2,3 an SU(2) index. (Anti)self- 
duality then means 0;*^^^°* = ±| jkUJ^^^^^ . 

Consider now a 2-form with values in the adjoint representation, -B^^. For this object we can 
consider the self-dual projections in both spacetime and algebra indices. Plebahski [8, 10] found 
that the two notions of self-duality coincide if the 2-form is constrained to satisfy the following 
quadratic condition, 

A B^ - -5'^B^ ABk = 0. (4) 
3 

This equation is called metricity constraint, and it is solved by B^ oc S*(e), where 

j:\e) = e'^ Ae' + ^e'jke^ Ae^ (5) 

is (twice) the self-dual part of the two form e A e. It satisfies^ 

J:' AJ:^ = 2ie6'U^x, (6) 

with e = e'^'^'"^e^e^epe^ the determinant of the tetrad, related to the determinant of the metric 
by det 5^,^ = — e^. The proportionality coefficient between B and S can thus be fixed to 1 
requiring 

B' ABi = 6ied^x. (7) 

Therefore a B field solving (4) encodes a metric via the tetrad field in (5). In order for these 
to give a real Lorentzian metric, reality conditions on B have also to be imposed. These are [10] 

B' A B*^ = 0, Im {iB' A B,^ = 0, (8) 

which imply respectively that B*'^ is proportional to the antiself-dual form, and that the deter- 
minant (7) is real. The B field itself can well be complex. These additional reality conditions 
are the trade-off for using the simpler self-dual variables instead of the full Lorentz group. 

In the following, we will refer to a i? field satisfying (4), (7) and (8) and the reality conditions 
as "metric". 

To come to the dynamics, let us also introduce an SU(2) connection a;^ - as an independent 
field, not the spin connection - and its curvature 2-form F^j,(cl') = duj^ + ^e^jkUJ^ Auj^ . Including 
the constraint (4) with a symmetric and traceless Lagrange multiplier '^ij, the action 



8^G^ 



B' A F,(c^) + i (m/,,- - ^-K5,j ]B'AB^ 



(9) 



^Note that also the antiself-dual part of e A e solves (4). The two solutions are distinguished by the sign of the 
right hand side of (6). 
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is equivalent, once the reality conditions for the fields are imposed, to {i times) the Einstein- 
Hilbert action for general relativity with cosmological constant A, plus a total divergence [8, 9]. 
It is instructive to see the equivalence at the level of the equations of motion. Next to (4), the 
other field equations are 

dB' + ie'jkio^ AB'' = 0, (10) 
F'{u) + (^'j - ^A5^)bJ' = . (11) 

When B is metric, (10) becomes the first Cartan structure equation, which - upon requiring 
invertibility of - identifies io as the (self-dual part of the) spin connection uj{e). This in 
turns implies that F*(u;(e)) is the (self-dual part of the) Riemann tensor (via the second Cartan 
structure equation). Finally, the 18 equations (11) have the following double role: ten of them 
give the Einstein equations for the tetrad, while five identify the Lagrange multiplier with 
the (self-dual) Weyl part of the Riemann tensor. The remaining three equations are gauge. See 
[2, 10, 16] and references therein for more details. 

Before moving on, let us add a brief comment on Plebahski's theory. The metric solution 
(5) to the constraints is unique if i?* A i?j 7^ is required, a condition that translates into the 
invertibility of the tetrad. Degenerate solutions with B^ A Bi = also exist, and they have no 
analogue in general relativity. The existence of this other sector makes Plebahski's gravity more 
general than Einstein's theory. Let us also recall that the canonical analysis of (9) gives the 
Ashtekar formalism [17], which is at the roots of loop quantum gravity [18]. 



3 Modified Plebanski gravity and non-metricity 

The modification introduced by Krasnov is to promote the cosmological constant term in (9) to 
a scalar function of the Lagrange multiplier itself, say $(\I'). The action then reads 



S[B, u;, ^ = j B'a FiiLo) + ^ + Sij f (^)) B' A B^ 



(12) 



As before, i = 1,2,3 is an SU(2) index, is symmetric and traceless, and the fields are 
complex. The reality conditions are the same as before. For the moment, we leave ^ arbitrary. 

When <I> = is constant we obtain Plebahski's gravity with A = — 3<^o; a non-constant ^ 
is responsible for deviations from general relativity. To understand the nature of the deviations, 
let us look at the field equations. The uj and B variations yield 

dB' + ie'jkUJ^ A = , (13) 
F'iio) + + <^{^)d'j^B^ = . (14) 

The key difference with Plebahski's gravity comes from the variation with respect to the field 
^. In order to take this variation, notice that only two algebra scalars can be built out of a 
symmetric and traceless ^ij, zi = Tr^'^ = ^''j^^i and 2:2 = Tr^'^ = ^'^j^'-'fc^'^j. Therefore $ 
is truly a function of zi and Z2- Taking also into account the tracelessness of ^, the variation 
gives 

f = 2 ^ + 3(^^^ - 1 ^ . (15) 

\ 0"^ J OZl \ 3 / dZ2 
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Thus instead of the metricity constraints (4), we now get the following equation, 

B' A B^ - ^'^B'' ABk = - ABk. (16) 

When the right hand side is non-zero, (5) is no longer a solution, and the theory is in- 
equivalent to general relativity. Hence the quantity 6^/5'^ controls the departures from general 
relativity. In the previous Section we recalled that when the right hand side of (16) vanishes, 
the solution B^ = 5]*(e) encodes a metric through (5). We are now going to show how a 
metric is encoded in B^ also when the right hand side does not vanish. 



3.1 Extracting the metric from the modified constraints 

The equations (16) are the key to understand the nature of the modification to gravity. They 
are manifestly not anymore constraints for B, but rather five equations^ relating B and . For 
instance, they can be used to determine the five ^^■''s as functions of the S*'s, once a specific 
choice of $ is made. The ^ so obtained can be then plugged in the equations (14) to obtain the 
dynamics of B. 

There is however an alternative procedure, as emphasised by Krasnov [2], that allows us to 
describe the theory in more familiar metric terms. This is based on the fact that the notion of 
Hodge self-duality naturally defines a metric (up to a conformal factor), even before imposing 
(4). Hence, simply declaring B to be self-dual suffices to endow the theory with a metric, and 
one can derive its field equations from the action (12). 

To be more explicit, consider the following two symmetric tensors, g^y and /i*-^, that can be 
constructed out of B^ [7, 15, 19], 

VW\9l. = ^e^^^'e.,,Bi^B^^^Bi, (17) 

B' A B^ = -B'' a Bk h'^ , /i\ = 3 . (18) 
3 

Here is the determinant of the Urbantke metric gj^^ introduced in [19, 10]. The crucial fact 
[7, 15] is that whenever h^^ is invertible, -B* is self-dual with respect to the Urbantke metric (17) 
- or to any metric related to it by a conformal transformation, because of the invariance of 
Hence the notion of self-duality together with B^ A Bi ^ suffices to have a unique metric, up 
to conformal transformations; we do not need to impose any constraints. 

But if a metric is already present, or rather a conformal class of metrics, what is the role 
of the metricity constraints (4) in the standard theory? It is instructive at this point to count 
the variables. A generic triple of 2-forms B^ has 3 x 6 = 18 components, of which 3 can be 
gauged away using the SU(2) symmetry. The remaining 15, thanks to the above result and 
when B^ A B^ is non-degenerate, can be parametrised in terms of the two tensors (17) and (18): 
_B* = B^{g^j^, K"^). Then, the role of the five equations (16) is to fix the five components of h^^: 
to constants in the standard case, to functions of '^''■^ in the modified case. 

*The trace of (16) vanishes identically. 
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To see how the fixing goes, consider first the standard case. Plugging (18) into the metricity 
constraints (4), we get h^^ = 5^^ . This kihs the extra five components in the {g^,y, h^-') parametri- 
sation of S*, and we are left with only the metric g^^ as the variable of the theory. An explicit 
computation then shows that g^^, coincides with the metric ej^elrjjj encoded in 5]*(e) via the 
Plebahski solution (5). The only conformal ambiguity left is the rescaling by a constant, but this 
can be fixed as in (7), or reabsorbed in the definition of Gn and A. So the standard constraints 
guarantee the existence of a preferred metric in the conformal class. 

In the modified case when the right hand side of (16) is non-vanishing, we can still use these 
equations to eliminate h^^ , but rather than being fixed proportional to the identity, this time 
its components satisfy algebraic relations with (the derivatives of) the ^'^•''s. Thus after solving 

(16) we have a new parametrisation B^{g^^^ ^^^). Plugging the latter in the torsionless condition 
(13) and solving this equation^ we obtain a (\I'-dependent) spin connection u}{g^ ,'^^^). Finally, 
inserting the solutions B{g^-' and u!{g^,^^^) into (14) gives us the dynamics for gj^^ and 

However, the conformal ambiguity is not trivial anymore in the modified case, for we can also 
consider metrics which differ from (17) by a ^-dependent conformal factor. One such metric 
arises if we solve the modified metricity constraints in terms of the metric encoded in S(e), 
instead of (17). As remarked in [2], this alternative construction is more convenient for practical 
reasons, and it is the one that we use in this paper. First of all, we restrict to the B^ A Bi ^ 
sector, and assume h^^ to be invertible, and thus diagonalizable since symmetric. Consider the 
matrix D^j of normalised eigenvectors, which is such that h^j = D'^ ^5^ i{D~^y j . h^^ has five 
independent components, hence also the diagonalizing matrix D. Call them c„, n = 1 ... 5. 
Using D(cn) we can bring (18) in a form reminiscent of the metricity constraint (4). Defining 
B^ = (D^^YjB^ , equation (18) becomes 

B' A & = ^ B^ A Bk 6'^ . (19) 

This equation is solved by B^ = S*(e), and therefore 

B' = D',{cn)^^{e). (20) 

Inserting this ansatz for B^ into (16) gives us algebraic relations between Cn and (the deriva- 
tives of) from which we obtain D^j{cn{'^)) = D^j{^). This gives a parametrization 
B' = B'{e^^,¥^) which solves (16). 

This way of constructing a solution to the equations (16) has the advantage of a simple link 
to the usual Plebahski form (5). Using this B to evaluate the Urbantke metric (17) one obtains 

~ /(^) ^li^iviJi where / is a scalar function obtained from contractions of three matrices 
-D*j(^). As anticipated above, the metric e^^e'lrju contained in (5) and the Urbantke metric 

(17) now differ by a non-trivial conformal factor which depends on ^ . 

Which of these two should be taken as the physical metric? As it turns out, there is no 
physically distinguished metric in the modified theory, and only the conformal class of the 
metric is specified.^ Thus both e^e^ry/j and g^^^, as well as any other metric in the same class, 

""That this is possible in general is shown e.g. in [20, 21, 22, 7, 15]. 

®For a discussion on this, see the literature by Krasnov and in particular [26]. 
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are treated on equal footing. This conformal ambiguity distinguishes the modified theory from 
Plebaiiski's. 

The ambiguity is expected to be resolved by matter couplings, and indeed in [26] a mechanism 
to distinguish a physical metric was shown, based on the requirement that test particles follow 
geodesies. In this paper on the other hand we consider only the vacuum theory, or its coupling 
to electromagnetism, which being conformally invariant is insensitive to the ambiguity. Thus 
the ambiguity is fully present, and it becomes relevant when discussing the singularities of the 
theory, because a metric has to be singled out in order to do so. To that end, we will consider a 
specific metric in the Sections on black holes below. We anticipate here that this is the metric 
^li^ivu which arises naturally in our construction of the spherically symmetric solution. 

3.2 Reality conditions and degrees of freedom 

One would expect the procedure just described to expose the presence of new propagating 
degrees of freedom with respect to general relativity. Remarkably, the Hamiltonian analysis 
performed in [5] (see also [11]) shows that this is not the case: also the modified theory has only 
2 propagating degrees of freedom. The key for the absence of new physical degrees of freedom 
is in the fact that the field equations for obtained from (14) are purely algebraic [7, 15]. 
The property is related to the peculiar complex nature of the action, and to the fact that reality 
conditions have to be additionally imposed, by hand. It is ultimately this imposition that it is 
responsible for the absence of extra degrees of freedom, both in the Plebahski theory (9) and in 
the modified theory (12). 

In this perspective, let us add a remark on the reality constraints before continuing. After 
the revival of Plebahski's complex formalism by Capovilla, Dell, Jacobson and Mason, it was 
later realized by Reisenberger [23] (see also [24]) that one can use real fields. In order to do 
so, one gives up the self-duality of the fields and takes B^"^ , uj^"^ and ij/^-^^^ in the full Lorentz 
group. The action is identical to (9) with Lorentz indices, and it is often called S0(3,l) (or 
non-chiral) Plebahski action. The absence of reality constraints is a manifest advantage of this 
formulation. 

The same modification proposed by Krasnov in (12) can be applied to the non-chiral S0(3,l) 
action. In its simplest form <I> oc Tr^'^, this was introduced by Smolin [14]. However, although 
the chiral and non-chiral Plebahski actions lead to the same classical physics, this might not 
be true for the modified actions. Above we mentioned how crucial the reality conditions are to 
ensure that the modified theory (12) only has two propagating degrees of freedom. The non- 
chiral action on the other hand does not have the reality conditions, thus the modification might 
include new degrees of freedom. This was recently pointed out by Alexandrov and Krasnov in 
[25] , where a counting of the constraints in the Hamiltonian framework showed the presence of 
six additional degrees of freedom. The modified non-chiral action thus appears to be a rather 
different theory, and its relevance for gravity still to be explored. 

4 Maxwell action 

The coupling of matter to the standard Plebahski action (9) was studied in [10], where actions for 
scalar, spinor and vector fields were proposed starting from a reformulation of the conventional 
curved-space matter lagrangians in terms of the fundamental fields and . This conservative 
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approach guarantees that when the B field is metric, the standard minimal couphng of matter 
fields to the metric is recovered. 



Such a reformulation is particularly natural for gauge theories, whose action can be written 
directly in terms of the spacetime Hodge star * = (l/2e)€^'^ pa- Thanks to this fact, the gauge 
theory lagrangian proposed in [10] can be immediately coupled to the modified action (12), and 
we will restrict our attention to this case. In units 47reo = 1, the Maxwell action reads 

S{g^,,A^) = ^ J FA*F = -r^ J d^x^gg>^Pg-'^F^M)FpM)- (21) 

To express this action in terms of the B field, a simple choice is^ [10] 

Sm[^. A,B] = ^J [if, B' A F{A) - i B' A S^] . (22) 

This is a first order formulation of electromagnetism, where the vector potential A is flanked by 
an auxiliary field ip. This action, which can be straighforwardly generalized to the Yang-Mills 
non-abelian case, is reminiscent of the first order formulation of Yang- Mills theory in flat space, 
sometimes called field strength formulation, or BFYM theory [21, 28]. 

When B is metric, (22) is equivalent to (i times) the Maxwell action (21). To see this, note 
that the vanishing of the variation gives 



A 



F{A) - ipjB^ 



0. 



(23) 



Recall that self and antiself-dual forms make a basis in the space of 2-forms. When the realitity 
conditions (8) hold, a self-dual B implies that its complex conjugate B* is antiself-dual. Then, 
a generic Maxwell field can be decomposed as 



F = F^'^^ + = F^^'^ B^ + F^^^ B*^, 
and thus B AF = B A F^+\ Hence (23) implies 

Fl,t\A)=^jBj,a. 
Substituting this result back into (22) we obtain 

Sm[^,A,B] = ±-jd^xe^^P-Fl+\A)F^p^\A) = ^jFA{l-i^)F. 



(24) 



(25) 



(26) 



Finally, if B^ is metric we can use (1) and get explicitly 



Su[^,A,B] = — 
lovr 



d'^x 



gg^''g^''FpM)FpM) 



-e^'^P'^F,,{A)Fpa{A) 



(27) 



'^The actions for scalar and fermion fields proposed in [10] can not be coupled to (12) the way they are, because 
they carry additional constraints, which are compatible only with the standard Plebariski ones (4) and not with 
(16) [4]. Hence suitable modifications of these actions have to be studied before scalars or fermions can be coupled 
to the modified theory. 

^Alternatives for the coupling of Yang-Mills theory to the non-chiral SO(3,l) Plebahski action have been 
considered in [14, 27]. 
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As promised, the first term of (27) is {i times) the Maxwell lagrangian (21). The second term is 
a total divergence, and as such does not contribute to the equations of motion.^ Neglecting this 
total divergence (and the similar one arising in the gravity sector), (9) plus (22) are equivalent 
to 



tot 



levrG J ^ 167r 
in agreement with the standard conventions [29]. 



(28) 



This proves the validity of the action (22) for coupling Maxwell theory to Plebahski's grav- 
ity. The constraints (23) are compatible with the metricity constraints (4) as well as with the 
modified constraints (16), thus the action (22) can also be coupled to the modified theory. This 
is what we do in this paper. That is, we consider (22) plus (12) as the total action for modified 
gravity coupled to electromagnetism. 

Since the Maxwell action (22) does not depend on the spin connection tu nor the field ^, 
the coupled action leads to the field equations (13) and (16) unchanged. The only equations to 
acquire a source are (14), which now read 

F'{lo) + {¥j + ^{'i^)6'j)B^ = SttGn t' (29) 

where we conveniently defined 

'^'S'matter _ ^ nupai (nr\\ 

6Bi ~ 4 ^ ^ 

pv ^ 

The algebra- valued 2-form r* = ^r^^dx^ f\dx'^ plays the role of the energy-momentum tensor in 
the standard metric formalism, and we shall refer to it as such. Notice however that unlike the 
standard energy-momentum tensor, r* is antisymmetric in the spacetime indices. For ^matter = 
given by (22) we get 

ri^ = -^^'Fl,-\A). (31) 
This satisfies ATi = 0, the analogue of the standard tracelessness condition. 



5 Spherically symmetric spacetime 

In the rest of the paper we consider the spherically symmetric reduced sector of the theory. 
This sector has been investigated by Krasnov and Shtanov [3, 6], who obtained the reduced field 
equations in vacuum. In this Section we review their results and extend them to the coupling 
to electromagnetism. 



5.1 Preliminaries: changing SU(2) basis 

Let us look back at the solution (5) of the standard Plebahski theory, which encodes a metric 
ds"^ = e*^ (8> e*^ — (8> — (8) — (8) via the tetrad e^. With the prospect of spherical 

^The second term is the standard CP-violating ^QCD-term, but with an i factor that makes it purely imaginary. 
As explained in [10], this complex action for Maxwell theory can be obtained from the usual one in the Hamiltonian 
framework through a complex canonical transformation, where the momentum conjugated to the connection is 
E + iB. This is analogue to what happens in the gravitational sector with Ashtekar's self-dual variables, where 
the momentum conjugated to the densitized triad is F -I- iK [17] . 
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symmetry in mind, it is convenient to work with a Newman-Penrose tetrad £^ = (/^, n^, m^, m^), 
which defines the metric in the form 

ds^ = I ® n + n ® I — m m — fh (^^ m, lf\nf\mf\fh = —ie d'^x . (32) 

The transformation between tetrad basis is / = {e^—e^)/V^, n = (e''+e^)/\/2, m = {e^+ie^) / ^/2, 
and 771 = m*. 

With the Newman-Penrose tetrad, the Plebahski solution (5) reads 

T} = I f\n — m f\rh , Yj'^ = n f\fh + I f\m , = i(n A m — / A m) . (33) 

This triple can be simplified with a change of basis in the internal SU(2) space. With respect 
to the standard basis Tj given by (one half times) the Pauh matrices, we consider the spherical 
basis To = rs, f± = ±(ri ±ir2). In this new basis we have, after a trivial permutation 123 i— > 312 
of the algebra indices, 

Y^ = lAn-mAfh, T,~^ =n Afh, =m Al . (34) 

A straightforward computation checks that this triple is still self-dual with respect to (32), and 
that the reality conditions (8) are satisfied provided / and n are real and rh = m* . Just as (33), 
the triple (34) provides a basis in the space of self-dual 2-forms. 

In the spherical basis fj, the Killing-Cartan metric is not diagonal any longer, instead 



2 

Accordingly we have 



T^irdj) = - , <5i, = I -2 I . (35) 

0-2 



Tr(S AS) = AS°-2S+ AS". (36) 
Consistency with (6) can be straightforwardly verified, 

6ij A S^' = (5,j A S^' = -6 / A n A m A m = 6 i e d'^x . (37) 

The non-trivial scalar product given by Sij is the trade-off for the simpler form of (34) . How- 
ever the advantage of the spherical basis fj is not only to simplify the triple of self-dual 2-forms 
from (33) to (34), but also to allow us to deal with tensor product of irreps using the familiar 
Clebsch-Gordan decomposition. This will be needed to describe ^ and the modified constraint 
equation (16). Before doing so, notice that the fj are a basis in the adjoint representation 1, 
but they are not normalised with respect to the scalar product (35). The transformation to the 
standard normalised basis |j, m) is: |1,0) = |fo), |1,±1) = |t±)/\/2. These extra factors of y/2 
have to be taken into account when writing the Clebsch-Gordan coefficients.^'^ 



'^"Alternatively, one can work with the normalised SU(2) basis. This, however, introduces factors of V2 in (34), 
complicating the field equations. 
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5.2 Solving the modified metricity constraints 

In order to solve the equations (16), we need to write explicitly the right hand side. To do so, 
recah that ^, which is a symmetric and traceless tensor product of two SU(2) adjoints, hves in 
the five dimensional irreducible representation 2. To find a basis in this space, we consider the 
(reducible) tensor product of two adjoints, 1(8'1 = 2©1©0. Here the irrep 2 is the symmetric 
and traceless part of the product, the 1 is the antisymmetric part, and the singlet is the trace. 
A basis in the irrep 2 can be obtained using the following Clebsch-Gordan decomposition, 

|2,±2) = |1,±1) ® |1,±1), (38a) 
|2,±1) = -^(|1,±1) |1,0) + |1,0) ® |1,±1)), (38b) 

|2,0) = -^(|1,1)©|1,-1) + |1,-1) ©|1,1)+2|1,0)©|1,0)) . (38c) 

In the case of spherical symmetry, the only non-zero component of ^ is along |2,0). Recalling 
the normalisation difference with the basis given by fj, we can write 

^' = -/3(f+ ® f_ +f_ ® f+ + 4fo fo) • (39) 

The minus sign is to have the same (3 as in [2, 3], where different conventions were used. 

To write the equation (16), we also need the identity in the spherical basis. The latter is the 
inverse 6^^ of (35), 

1 = (5*^fi (g) fj = fo To - -f+ f_ - -f_ (g) f+ . (40) 

This identity naturally belongs to the singlet and corresponds to the standard SU(2) scalar 
product g"^" = {—l)'^6m-n, once the different normalisation of the f± is taken into account. 

We now have all the ingredients to write the metricity equations (16). Using the metric 6ij 
in (35) to take the scalar products, we compute 

¥k'i>''^ fi(g>fj = /3^(-2f+(g)f_ -2f„ 0f+ + 16fo(g)fo), (41) 
^«^$fc,^'if. 0f^. = -/33(4f+0f_+4f_(g)f+ + 64fo0fo), (42) 
from which zi = ^'^^^'-'j = 24/3^, Z2 = ^-^-^^i^'j = -48/3^. 

Being both zi and Z2 functions of /3, it follows that the modification to general relativity 
is characterised in the spherically symmetric sector by the functional ^{P) of a single scalar 
function. As a consequence we also have that the partial derivatives are simply 6zi^ = $/3/48/3, 
Sz2^ = — *l'/3/144/3^, where we have introduced the shorthand notation <I>^ = 6^/6(3. Notice 
that $ and ^ have the same dimensions, thus is dimensionless. 

Putting all this together, eq. (16) gives 

i^^ - ^) (^+ ® ^- + r_ f+) + ($/3 + 1) fo © fo 

Prom this we can obtain the following five independent equations, 

B+ AB+ = B- AB- = 0, (44a) 
A S° = 5" A S° = , (44b) 

2B+ AB' +B^ AB^ = ^$/3 B^ ABk. (44c) 



B' A B^ fi ®f.j=\B^ A Bk 
3 



(43) 
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Using the scalar product (35), the last equation can be written as 



= ^/^-^ (44d) 
BO A BO 4(1 + $^) ^ ^ 

Proceeding as discussed above in Section 3.1, we now look for a solution in the form 

B' = D',{cn)^^e), (45) 

with S*(e) given by (34), and c„ five parameters to be related to through (16). In the 
spherical basis, a generic linear combination that preserves self-duality is given by (see e.g. 
Appendix of [2]) 

/ CO \ 

D{cn) = c+ c+_ . (46) 
V c_+ c_ / 

Plugging this into (44a), we have 

n n C+C-+C+-C-+ - 2 

= ^-^^+ = 0' -2(co)^ = 4(rT¥^- 

A simple parametrisation of the solution is given by c_| = c ^ =0, Cq = 1 and c+ = c— = c 

with 

(48) 



2_ 



2(1 + ^/3)' 

In conclusions, the i?* field given by the triple of 2-forms 

B^ = cn/\fh, B^ = cmAl, B^ = lAn — mAm, (49) 

solves the metricity equations (16) and the reality conditions, and it is self-dual with respect 
to the metric (32). Note that although might be negative, hence B complex, (32) is always 
real and Lorentzian (up to singularities). Note also that the relation between B^ and the metric 
(32) breaks down when <I>^ = 2 or = —1. These situations can indeed arise, and are the 
non-metric singularities which we will discuss below. 

Needless to say, having a metric at disposal is instrumental for the physical interpretation 
of the theory. The strategy is to use (49) to write the remaining field equations as differential 
equations for /, g and /?. Once the theory has been reformulated in such more familiar terms, 
also the effect of the modification become more transparent. 

Thus far, the only place where we used spherical symmetry was the fact that ^ has the 
single component (39). (In the general case the equations (44a) have non- vanishing right hand 
sides and depend on the whole five components of ^, which get related to the coefficients c„). 
We now specialize the Newman-Penrose tetrad, and thus the solution (49) for the B field, to 
the case of spherical symmetry. We choose spacetime coordinates as follows: 

11/? 

I = — ^ {f dt — g dr) , n = — ^ [f dt + g dr) , m = {dO + i sin 9 dcp) , (50) 
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and m 



m 



Through eq. (32), they define the metric 



fdt^-g'^dr^-R^ide'^ + sm^edcj)'^), e = fgR^s\n9, (51) 

which can be recognised as the standard ansatz for a spherically symmetric spacetime. We 
initially take f{t,r), g{t,r) and R{t,r) functions of time as well as the radius. 

Using the tetrad (50), we can write the components in (49) in the following matricial 
form, 

\ / ±Rf -iRf sine \ 



B' 



/ /<7 



iR^ sine 
\ 



B 







Rg 




^iRg sin e 





(52) 



5.3 Maxwell fields 

With the explicit form of B^ at hand, we can immediately compute the Maxwell stress-energy 
tensor (31). It is straightforward to show that when a Maxwell field respects the spherical 
symmetry, the field strength takes the form 



F{A) = E(t, r) / A n + B(t, r)mAm 



(53) 



with two scalars, E{t, r) and B{t, r), corresponding respectively to the electric and magnetic field. 
Then, it follows from the source-free Maxwell equations dF = and d* F = (which are still 
valid in the modified theory in terms of the derived metric) that B[t,r) = const, x R[t,r)~'^, 
E{t,r) = const, x R(t,r)~'^. For simplicity, hereafter we restrict our attention to the case 
B(t,r) = 0, and thus in units Aireo = 1, 

Q 



F{A) 



I An. 



R{t,ry 

where the constant Q is chosen to be the conserved electric charge defined by 



(54) 



(55) 



on a 2-sphere 5^. 

To evaluate eq. (31) we need the self and antiself-dual parts of (54), which can be computed 
from (24) and the explicit form (52) of B. The only non-vanishing components turn out to be 



F, 



(±) _ Qfg 

2i?2 ' 



f!;^^ =±i^sin0. (56) 

(57) 

Finally, plugging (56) and (58) into (31), the only non-zero components of the energy-momentum 
2-form compute to 

Q^fg . 



tr 



±1 — sm ( 
2 



2ip+ B^.-lip- S 



From the equation of motion (25) we have 
which gives 



Q 



8ttR^' 



'8ttR^ 



sm ( 



(59) 



15 



5.4 Cartan equations 

The next step is to solve the Cartan equation (13), using the exphcit form of in (52), to 
compute the connection to in terms of the metric and (3. Since we are using the spherical basis, 
we have to adapt the covariant derivative. Recalhng its general definition d^B = duj + [^,B] 
and using the commutators in the fj basis, we get 

dB° + 2(0;- AB+ -UJ+ AB-) =0, dB^ ± (uj'^ A B^ - uj^ A B°) = . (60) 



To solve the equations we proceed as follows. First, we use (52) and write down the non-zero 
components of the 3-forms dB^. Then, we take the wedge product of the equations (60) with all 
four basis vectors Z, . . . ,m, and read off the components of to in this basis, as functions of the 
metric and /3. The procedure is straightforward but space consuming, so we report here only 
the final result (more details can be found in [2, 3]), 



± 



1 



P±{d9 =F i sin( 



F dt + G dr - icost 



(61) 



where we have introduced the following shorthand notation. 



9* J* 



F 



{Rf*y R'f* 

Rg* 



c^Rg* ' 



G 



{Rg*y Rg* 



Rf* c^Rf. 



(62) 



with /* = c/ and g*, = eg. Here and in the following ' and ' denote, respectively, the derivative 
by t and r. 



The curvature in the spherical basis is 
F° = (ia;° + 2a;~ A u;+ , 
where from (61) we have 



( F G ^ \ 

F' G' 

isg 

V y 

/ 

P+P- 

4 1 

V ISO 



Here sq = sin0 and cq = cos 9. 






-is0 

o2 



F^ = du^ ±uj° Auj^, 



± 



duuj 



/ -P± ±i se P± \ 

-P4 ±iseP4 

±ic0P± 
V J 

/ -F ±iseF \ 

P± -G ±is0G 
2 

\ ice ±C0So I 



(63) 



(64a) 



(64b) 



5.5 Field equations 

Using (39) and (63), the field equations (14) in the spherical basis read 



dw^ ± A + (/? + $) B^ 
duP + 2w~ A CU+ + ($ - 2/3) P° 



SttGn ■ 



(65a) 
(65b) 
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Here co = uj{B,f3) through (61) and $ = as seen above. For a spherically symmetric 

configuration of Maxwell's field, we found that the only non-zero components of the energy- 
momentum tensor are r^J. and r^^. Using also the explicit expressions (52) and (64), the non- 
vanishing components of equations (65a) are 



P±±P±Ft{P + $)i?/, = , (66a) 
P4 ± P± G - (/? + $)i?ff, = , (66b) 



and those of (65b) 



G-F' + i^-2p)^ = SttGn T^r , (66c) 

1-P+P^+R'^{<^-2P) = 8ttG^-^tII^. (66d) 

I sm d ^ 

Among the eight equations (66), only three are independent because of the Bianchi identities 
[3]. Below in Section 5.7 we write the three independent equations in the form suggested in 
[3]. Before doing so, we show the validity of part of the Birkhoff theorem that any spherically 
symmetric vacuum solutions are static. 

5.6 Birkhoff's Theorem 

In [3], Krasnov and Shtanov proved the Birkhoff's theorem that the spherically symmetric 
solutions of this theory are necessarily static for the vacuum case. Their proof is slightly un- 
conventional, so we find it useful to reproduce the same results for our case, following a more 
standard procedure [30]. We assume that R{t,r) be not constant, dR / 0. This, in particular, 
implies that we can choose the gauge so that R{r, t) = r. Then, the field equations above reduce 

to 

P± = -, F=^ + A('i_^y G=f^. (67) 
9* 9* rg^ \ J 

The first two of equations (66) simplify to 



\g,±-F^{fi + ^)rU = Q, (68a) 

gi 9* 

— 9*^-r- 

gi J* 9' 



-g'^±-Ml-G-{(3 + ^)rg, = 0. (68b) 



We now show that the solution is necessarily static, i.e. that f = g = c = 0. First of all, 
summing the two equations in (68a) gives g^, = 0, thus G = from the third of (67). It then 
follows from (68b) that /? -1- $ is a function of only r and thus (1 + ^is)^ = 0. As we briefly 
mentioned above <I>/3 = — 1 corresponds to the singular case, which we are not going to consider 
here. We therefore conclude that (3 itself is t- independent and so are $ and c^. Subtracting the 
two equations (68a) gives -F//* = rg^,{(3 + $), implying that {F/ f^,)' = 0. Then, it follows from 
(67) and {F/ f^)' = that (/^//*)" = 0, namely /^//* = a(r) with a{r) being some function of 
r; the latter equation allows for a time-dependent integration constant, f^:{t,r) = b{t) J a{r)dr, 
but this is of the type that can always be absorbed by a change of coordinates t i— > r with 
dr = h{t)dt. 
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Note that besides this staticity property, Birkhoff 's theorem also asserts that such a spheri- 
cahy symmetric static metric, if satisfying the vacuum Einstein equations, is necessarily locally 
isometric to the Schwarzschild metric [30]. For our case, in order to have the corresponding 
assertion, we have to integrate eqs. (66c) and (66d) with the stress-tensor for non-vanishing 
Maxwell fields. We perform this integration only numerically and do not have analytic solu- 
tions. For this reason we do not have the corresponding assertion in the same decisive level as 
the original Birkhoff 's theorem for the vacuum Einstein gravity. 

5.7 Equations of motion for static fields 

Using the staticity and the choice of gauge R{r) = r, the equations reduce to 



rf*9* rgl 



P + ^ = ——F = (69a) 



- F' + ($ - 2/3) ^ = SttGn t° (69b) 
1 - ^ + r2 ($ - 2/3) = SttGn rO (69c) 

Here (69c) is an algebraic relation, and (69a), (69b) are differential equations for /^,, g^, and 
f3, of which only two are independent thanks to the Bianchi identities. Following [3], we now 
reexpress these equations in a more convenient form. First of all, (69a) and the second of (67) 
immediately imply that 

%^ = -ifi-i). (70) 

j*g* r \ & ) 

Secondly, taking the derivative of (69b) and using the above equation, we obtain 

/3'(2 - $^) = -- (3/3 + ^ SttGn tI) . (71) 

The three equations (69c), (70) and (71) are the field equations of the theory in the spherically 
symmetric sector. They give dynamics for /=„, g^, and (3, and thus for the metric (51). For 
(5 = these equations reduce to Eq. (35) of [3]. Finally, using the explicit expressions of the 
energy-momentum tensor (59) and of (eq. (48)), and defining x(r) = f*{r)g^{r), we have the 
independent equations 



^ - f72a) 
X - r(2-c|.,)' (^2^) 

{2-^,)P' = _^ + ^, (72b) 
1 = l-i2P-^y-^. (72c) 

Before proceeding, let us check that for constant ^ = —A/3 we recover the standard results 
of the general relativity case. 
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5.8 Recovering the standard Reissner-Nordstr0m metric 

When <I>^ = we have = 1 and, from (72a), that f^^g^ = fg is a constant. Taking the latter 
to be 1, / = 1/g. We are left with a single differential equation, from (72b), 

^' = -^ + ^. (73) 



This equation is solved by 



Pin = — 3 — , (74) 



where M is the integration constant. Finally from (72c) we read 



— = 1 + 2 ('^5) 

gZ ^ ^ 



which can be immediately recognised to give the 500 = 5 ^ component of the Reissner-Nordstr0m 
metric with cosmological constant A. 



5.9 Metric singularities 

The field equations (72) give the dynamics for f^, g^ and (3, in terms of which the metric (51) 
reads 

ds' = [f^ dt' - gl dr^] - d'n . (76) 

Notice the possibility of singularities if values of r exist such that = 2, —1. These are the 
non-metric singularities found and discussed in [3]. As shown above in Section 5.2, these are 
points where the relation between B and the metric breaks down. To understand the theory, it 
is useful to distinguish the following three types of apparent singularities which can occur: 

• Points where gj(r) diverges and x{r) is finite are Killing horizons, the outermost of which 
may be viewed as the black hole (event) horizon, depending upon the asymptotic structure. 

• Points where = 2. We call this non-metric singularities of type (1). 

• Points where = —1. We call this non- metric singularities of type (2). 

The nature of these non-metric singularities will be discussed below, with explicit solutions at 
hand. 

It is at this point important to comment on the conformal ambiguity present in the theory. As 
discussed earlier in Section 3.1, the procedure to extract a metric is defined up to a conformal 
transformation, and in the electrovacuum theory there is no principle which would select a 
physical metric in the conformal class. In particular, multiplying (76) by a conformal factor 
generically depending on (5 gives an equally valid line element in this modified theory of gravity. 
A possible mechanism to distinguish a preferred metric has been studied in [26] by introducing 
test particles and demanding that they follow geodesies. 

Since we restrict our analysis to the coupling to electromagnetism and Maxwell's equations 
are conformally invariant, the ambiguity is present in our case, and we have to arbitrarily select 
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a metric in order to investigate its singularities. Thus in the absence of a physical criterium to 
distinguish a preferred metric, we choose to work with the the line element (76), which emerged 
naturally in our construction.^^ 



6 Choices of $ 

Thus far we have kept $ arbitrary, without worrying about the origin of its non-constancy. 
Different choices of $ have an impact on the classical behaviour of gravity in different ways and 
at different scales, allowing one to entertain the possibility of explaining in this way current 
puzzles, such as anomalous rotational curve in spiral galaxies, without appealing to dark matter 
scenarios. To that end, tailored profiles of $ have been investigated in [3, 6]. 

In introducing this modification of gravity, Krasnov pointed out that it could be seen as an 
effective action coming from a model of quantum gravity where the fundamental field is not the 
metric, but rather the 2-form B. Recall that $ has length dimensions = [^'] = —2, and that 
it is a function of zi = Tr ^'^ and Z2 = Tr with dimensions [2:1] = —4 and [Z2] = —6, we 
see that a non-trivial dependence of <I> on ^' introduces coupling constants with non-zero length 
dimensions. A natural coupling constant with such dimensions is the Planck length, which leads 
us to postulate a quantum mechanical origin of a non-constant <I>. At present, it is not clear 
what quantum gravity model could lead to such an effective action, but the use of B instead of 
the metric as the fundamental field suggests to look at spin foam models as possible candidates. 

In this paper we do not attempt to quantize the theory, we are only interested in classical 
consequences of different choices of ^. Consider for instance the case in which $ is 'analytic' 
near the origin in the sense that it can be Taylor expanded as 

$(^') = $0 + aiTr^'2-ha2Tr^'3-F... (77) 

As [^] = —2, it follows that [an] = 2n, e.g = 5„ (?iGn)" with 6„ dimensionless. Lack of control 
on the possible origin and on the quantum properties of the theory makes it hard to restrict the 
admissible forms of $, and a priori we have no reasons to demand analyticity in the sense of 
(77). On the other hand, by studying the implications of a non-constant ^ at the classical level 
we can hopefully obtain restrictions on the admissible forms of 



6.1 On asymptotic flatness 

A useful requirement is that in the modified theory, there exist such solutions that approach at 
large distances an asymptotically flat solution of general relativity. This can put restrictions on 
the admissible choices of Taking A = for simplicity, we demand that 

c2 ^ 1, (78) 
f, ^ l + 0(l/r), (79) 
f3 ^ 0, (80) 

in the limit of large area radius r — > 00. Since = (2 — <I>^)/2(1 -|- ^p), the first of the asymp- 
totic flatness condition, eq. (78), implies that cannot approach a non- vanishing constant, and 

^^Another natural choice considered in the literature [3] is to work with the metric such that the volume form 
coincides with B'^ A Bk/3. This amounts to multiplying (76) by (1 + ^p)'^^^. 
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hence that with eq. (80), $ cannot be Unear to /?. Therefore the asymptotic flatness conditions 
above require that if Taylor-expanded in P, $ must start from at least the quadratic order of /?: 
$ oc hG^0^ + 0(/?^). This also is a justification of our discussion about the form, eq. (77), from 
a physical view point. The simplest modification to general relativity is thus to truncate the 
expansion (77) to the first non-trivial term, $ = $o + oi Tr = + 24ai/3^. (When $o / 0, 
the solution would be asymptotically (anti-)de Sitter.) This quadratic case will be discussed in 
detail in the subsequent sections. 

Since our present metric is static, admitting the time-like Killing vector field = {d/dty\ 
the ADM mass may be given by the Komar integral on 2-sphere, 5^, at spatial infinity, 



M = -^[ *(Ve). (81) 



(Note that the conserved electric charge, Q, is always given by (55), irrespective to the conditions, 
eq. (79).) Since 



f 2 _ — 2 



1 - (2/3 - -/32)r2 + 0(l/r2) 



(82) 



x2^1 + 0(l/r), p^- + 0{l/r^), (83) 



the second of eq. (79) implies that 

C 
r 

with some constant C. (In fact, integrating eq. (72a) at asymptotic region implies ~ 
exp[-(AC)/r^] ~ l+0(l/r^).) Then, the ADM mass — which is now written M = r'^c^{p)' /2x — 
is well-defined and becomes, C = G^M, in accordance with the standard GR results, eq. (74). 

6.2 Linear modification 

When $ is linear in /?, then x would either diverge or vanish at r — > oo. Therefore solutions 
in a modified theory linear in (3 would not be asymptotically flat. Since linearity in (3 can only 
happen if ^{^) is non-analytic at the origin in the sense of the expansion eq. (77), one may ask 
whether the lack of asymptotically fiat static solutions characterises all "non-analytic" choices 
of <I>. This would be a rather useful restriction on the admissible choices of <I>. We leave this 
question open for future investigations. 

Nevertheless, the linear modified case can be solved easily, and thus it is instructive to 
consider this case first. Let us take $ = $0 + with a being some dimensionless parameter. 
In this case is constant and eq. (48) reads 

2 ^ (84) 
2(1 + a) ^ ' 

The constancy of (? trivialises some of the non-linearities of the field equations, thus significantly 
simplifying the study of the spacetime structure. Assuming a 7^ 2,-1, the field equations are 

x' 3a 



(2 - a) r ' 



(85a) 



(2-a)/3' = -^ + ^, (85b) 



l + <^y-{2-a)l3r^-^^. (85c) 



gl 
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The equation (85a) can be easily solved to give 

x{r) = ci = ci r^^-^')/'^' , (86) 
with ci being an integration constant. From (85b) we obtain 

/3(r) = C2r-^ + ^g^r-^ (87) 

where C2 is another integration constant. The integration constants can be fixed to reproduce 
the standard GR result (74) for a = 0, specifically ci = 1, C2 = G^M. Then we have 

r^i^ 2a — 1 

Note the different asymptotic behaviour from the standard Reissner-Nordstr0m solution, and 
in particular the absence of asymptotic flatness. This means, in particular, that M cannot be 
interpreted as the ADM mass. Note how the location of the zeroes of depends also on 
a. In the next Sections, when we treat the quadratic modification, we will see that a similar 
modification of the horizon structure happens within solutions which are asymptotically the 
Reissner-Nordstr0m metric. 

7 Modified self-dual gravity with Maxwell field: Quadratic case 

In the rest of the paper we focus on the quadratic modification $(/3) = $0 + 24ai/3^. For 
commodity, we redefine ai = hG^A/48, and from now on we work in units h = Gn = 1. We 
then have $(/3) = $0 + A(3'^/'^, = Aj3 and 

2 -A/3 



2(1 + ^/3) ' 

which is positive for A(3 S (—1,2) and negative otherwise. 

As a consequence of this choice of <I>, we expect larger departures from general relativity in 
regions of larger curvature. 

7.1 Nariai-Bertotti-Robinson solutions 

Our general formulation of the field equations (66) allows us to look for solutions with R{t, r) = 
Rq constant. Such solutions exist in general relativity: They are not asymptotically flat and have 
a Cartesian product structure. The simplest case is the neutral Nariai metric [31], a solution to 
the Einstein equation with a positive cosmological constant given by a direct product d«S'2 x S^. 
Setting R{r,t) = Rq in (66), we have 

P± = 0, F=^, G = ^. (90) 

g* J* 

Equations (66a) and (66b) reduce to 

(/? + $(/?))/, = (/3 + $(/3))5* = 0, (91) 
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whereas (66d) reduces to 

l + Rl{^{[5)-2[5)=^. (92) 

Avoiding the unphysical case of /=„ and everywhere vanishing, equations (91) give an algebraic 
constraint for (3. Thus, the solutions to the constraint (3 + ^(/3) = will have constant Weyl 
curvature. 

In general relativity, i.e., $ = $o = — ^/3, the constraint is 13 = A/3, which plugged into 
eq. (92) gives 

Rl = \. (93) 
This fixes the value of Rq dynamically, and the Nariai solution follows from (66c). 
In the quadratically modified case, the constraint is 

i^/32 + /?-iA = 0, (94) 



from which (3 = {-l± yjl + {2/2,)KA)/A. Substituting this value into (92) we have 

Rl = , = . (95) 

±Vl + (2/3)AA - 1 

Therefore the Nariai solution still exists, but with a different value of its constant radius. This 
new value now depends on the modification parameter as well as the cosmological constant 
A. For A > 0, we take the solution of /3 with positive sign, and iig (95) is positive provided 
A is positive, as in (93) above. For A < 0, the solution with positive sign again requires A > 0, 
but the solution with negative sign allows any A < 3/21^41 . This, in particular, implies that the 
Nariai metric of dS^ x 5*^ can be a solution even when A < if A < with negative sign of 
eq. (95). It is however only the solutions with positive sign and A > that reduces to eq. (93) 
in the limit A 0, i.e., the limit to the standard general relativity case. 

Another case is the Bertotti- Robinson metric [32, 33], a solution to the Einstein-Maxwell 
equations with vanishing cosmological constant, given by the direct product of AdS2 x 5^ with 
each 2-dimensional manifold having the same constant curvature radius. This metric has van- 
ishing Weyl curvature and therefore is still a solution to the modified theory. Indeed, taking 
A = in (94) we see that /3 = is a solution also for A ^ 0. 

7.2 Perturbative analysis of charged black holes 

Let us now look for asymptotically flat solutions. Fixing R{r) = r, the equations (72) take the 
form 



x' 1 3A(3 



X r 2-AI3' 



(96a) 



P'i2-Ap) = -^ + ^, (96b) 



J_ = i-r^(^2P-lA(3^-^o)-%- (96c) 



With respect to general relativity, two immediate effects of the modification are visible: 
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• Prom (96a), x =const is not anymore a solution, so we do not have the famihar form, 

= 1/5^, of the sphericahy symmetric Einstein-Maxwell system. 

• In eq. (96b) the term proportional to A makes the differential equation non-linear. 

The non-linearity of the system is still manageable for Q = 0, the case that was solved exactly 
in [3], and which we review below. However for Q the situation is more involved, and 
unfortunately the system eluded our attempts to find exact analytic solutions. A solution can 
be found perturbatively in the regime \Ap\ <C 1, giving 

«0 ^ --^.^(m4)%„(.^, 

where we fixed the integration constants to match general relativity at zeroth order in ^. A 
posteriori we see that the approximation is valid for r ^ (AM)^^^. The perturbative solution is 
useful to fix the integration constants, because it turns out that, unlike the linear modification (cf. 
(87)), the quadratic one maintains the right asymptotic behaviour. This important property was 
obtained for the vacuum theory in [3] , and here we extend it to the coupling to electromagnetism 
by a perturbative analysis. Below in Section (9) we show it also using the non-perturbative 
numerical solution. Before studying numerically the system, let us discuss the vacuum case 
Q = 0, that can be solved analytically as found in [3]. 

8 Vacuum black holes revisited 

In [3], Krasnov and Shtanov solved the system (96) for Q = and found an analogue of the 
Schwarzschild black hole. Here we review their results and extend them to include the locations 
of the event horizon. The knowledge of the analytic solution allows us to test our numerical 
integration code, which we use for the plots of this Section and the next Section on the charged 
case. 

The Krasnov-Shtanov solution is given by 

x{r) = e-'^^f^ , (99) 

with /3(r) implicitly defined by 

r(/?)=(f)%^^ (100) 

M being an integration constant, corresponding to M in eq. (98).^^ From (96c) with Q = we 
^■^This solution can be expressed in terms of the Lambert W function (see e.g. [34]), 

«'-'-l-(-^)' 

We thank Robert Mann for pointing this out to us. 
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(97) 
(98) 



can immediately extract the non-trivial components of the metric, 



1 ^2 2(1+^/3) 



5oo(r) = ^f't= ' ^"^^ 1 - (2/? - ^Ap' - '^o)r' , (101) 



1 2 2(1 + ^/3) 



1 



2 



c2 ^* 2 - 



1 



9u{r) = -^9t = - ~\^ ' 7a' 1 - (2/3 - - $o)r" . (102) 



2 



One can easily check that setting ^4 = gives back the standard Schwarzschild-de Sitter metric 
representing a black hole of mass M. 

At values of r where 1 — (2/3(r) — ^A/3^(r) — ^Q)r^ = 0, we have a coordinate singularity that 
can be removed via a change of coordinates. If > outside, then the coordinate singularity 
may be interpreted as the event horizon. The metric becomes degenerate also at the points such 
that $/3 = A(5 = 2,-1. These are non- metric singularities where — > (type (1)) or — > 00 
(type (2)), and where the relation between the fundamental field B and the metric breaks down. 
To see whether and when these singularities occur, note first that the solution for /3(r) obtained 
inverting eq. (100) is always positive for r > 0. This means that A(3 = 2 can only occur if ^4 > 0, 
and the one at Af3 = —1 if yl < 0. Recalling that A = — 3<I>0) ^ positive (resp. negative) A has 
a coherent sign with a negative (resp. positive) cosmological constant. 

Note that $0 enters the equations in the same unmodified way that the cosmological constant 
does in general relativity. In the following we will restrict attention to <I>o = and focus on the 
effect of the modification introduced by a non-constant <I>^. 

What is the nature of the non-metric singularity? At these locations, the metric becomes 
singular or divergent, see eqs. (101-102). However it was shown in [3] that these are not physical 
singularities for the theory: a change of coordinates exists such that all the fundamental fields B 
and Lu are finite. To see supporting evidence that the singularity is not physical, we can also look 
at scalar invariants. One simple scalar invariant can be constructed out of the field equations 
(14) as follows: first, contract the equations with the wedge product of B^. This gives a scalar 
density. Second, divide by the volume element provided hj V = ^B^ A Bi. In the standard 
general relativity this procedure would give the Ricci scalar. In the modified theory, it defines a 
scalar quantity S = B^ A Fi/V . Taking $0 = for simplicity, the equations (14) give after some 
algebra 

S = \Af3\ (103) 

This shows explicitly the finiteness of this scalar quantity at the non-metric singularities of both 
type (1) Af3 = 2 and (2) A(3 = —1, and justifies, in part, the claim that these non-metric 
singularities are not physical singularities of the theory. What is happening is not that the 
theory is breaking down, but rather at the locations of these non-metric singularities, there 
doesn't exist any metric with respect to which the 2-forms, eq. (49), is self-dual. 

The two cases of positive and negative A are qualitatively different (see Fig. 1), and it is 
convenient to treat them separately. 



^''Note that here and in the following, by goo and pn we mean the {t, t)- and (r, r)-coordinate components of 
the metric; they are not the components with respect to frames, {e^}, or {l^}, introduced in Section 5. 
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Figure 1: Paradigmatic plots of r{(3) from (100), for A > (red) and ^ < (blue). In the first case a 
minimal possible radius arises. In the second case the function is always invertible. 



8.1 Positive A branch 

When A is positive we see from Fig. 1 that r{j3) in eq. (100) has a minimum. Consequently, 
there is a minimal value of r that can be reached. Imposing the derivative of eq. (100) to vanish, 
this minimum is found at Tnm = (e^M/2)^/'^. We call this point Vnm' because it is precisely 
where Af5 = 2, thus at this point there is a non- metric singularity of type (1). (3 can be inverted 
either in the domain r G (r^MjOo), or r G (0, Tnm). The first one has the standard asymptotic 
behaviour (5 ~ (r i— > oo) of general relativity, and we look at this case first. In this branch 
we have A[3 G (0, 2), thus (? is always positive; consequently, the outmost zero of fl corresponds 
to the event horizon, and changes in the signature of the metric come from changes in the sign 
of only. It is natural to ask whether the singularity at r^M hes within the event horizon, or it 
is naked. The horizon is located at /^(rn) = 0, but this equation cannot be solved analytically. 
To visualise the situation, we find numerically the zeroes r-fi{A) of and compared them with 
ri^yi{A). See Figure 2. To further illustrate the structure of the spacetime metric, we plot /3 and 
for various values of A in Fig. 3. The complete behaviour of the metric can be easily inferred 
from these plots and eqs. (101) and (102). 

We see from Fig. 2 that both rYi{A) and r^yi{A) increase with A. The radius of the non- 
metric singularity grows in A faster than the horizon radius, and there is a value of A for which 
the radius of the singularity reaches the horizon radius. This value can be computed by requiring 
f^{Al3 = 2) to be non-positive (see Fig. 3), and one finds A < 2e^M^. This means that for 
fixed A there is a minimal mass, Mq = y]4/(2e^ for which the black hole exists: for smaller 
values of the mass, the non-metric singularity is naked. 

At the non-metric singularity both qqq and gn are diverging, but this is not a physical 
singularity as we have discussed above. Hence the theory itself is not breaking down and the 
solution can be continued to describe a "spacetime" beyond this point. From the point of view 
of the metric, continuing the solution means using the other branch of the inverse of r(/3) in 
eq. (100). This is the branch where Ap G (2,oo) and we move from r = tnm to r = oo as 
we recede from the non-metric singularity. Since is always negative there, it follows from 
eqs. (101) and (102) that the behaviour of the metric inside the non-metric singularity is similar 
to the case of Kantowski-Sachs cosmological spacetime, with the singularity at r ^ oo (and thus 
(3 — oo, goo 0~, — > 0^) being a true curvature singularity. We do not report here plots 
of the solution in this branch. The overall picture can be summarised in the Penrose diagram 
already present in [3] (e.g., for the region of Kantowski-Sachs spacetime, see VI and VIII in Fig. 
1 of [3]), where further discussion of the global structure is also given. 
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Figure 2: Location of the horizons and non-metric singularities as functions of A, for M = 1. The red 
continuous curve is the non-metric singularity Tnm for A > and the green continuous curve for A < 0. 
The blue interpolated dots are numerically found zeroes of f1 representing the radius of the event 
horizon. For negative A, we have two horizons. The values of A which undress the singularity compute 
analytically to A = 2e^ ~ 14.78 and A = — 125/(8e) ~ —5.75. Analytically we can compute that the two 
zeroes of present for negative A merge at A = — 54e^^ ~ —7.31. Yet notice that the zeroes numerically 
found do not match exactly at this value: this shows explicitly the limitation of the numerical integration. 
The error can be estimated to be ~ 5%. 



A>0 




Figure 3: Plots of (continuous) and 2 — Af3 (dotted) for different values of ^ > 0, with M = 1. From 
left to right, A = 0.2 (brown), A = 1 (blue), the limiting case A = 2e^ (green), and A = 50 (red). The 
dashed black line is the standard general relativity result, drawn here for comparison of the asymptotic 
behaviour. The zeroes of the dashed lines are the locations of the non-metric singularities A/3 = 2, where 
= and this branch of the solution interrupts. 



8.2 Negative A branch 

When A is negative, the function r(/3) in eq. (100) is invertible for any r > 0. Since Aj3 £ 
(—00,0), the A/3 = 2 singularity cannot occur, but rather the A/3 = —1 one. At large radii, 
P approaches zero with the standard r~'^ behaviour. As we decrease r, we reach the value 
Tnm at which Ap = — 1 and the type (2) non-metric singularity occurs. The location is easily 
determined plugging A/3 = —1 in (100): we find r^M = {M\A\y/^e^^/^ . There, = 0, and 
both goo and gu vanish. As already seen, this is not a physical singularity of the theory. The 
solution continues after this singularity (recall /3' never vanishes, so eq. (100) can be inverted 
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for all r > 0), with < and the metric signature is switched to — I as (3 becomes large. 

Again, we can numerically find the zeroes of to see whether the type (2) non-metric 
singularity is hidden or not by the event horizon, see Fig. 2 for ^ < 0. In Fig. 4, we plot and 
p. The plots show that approaches the Schwarzschild solution at large distances, and hence 
we can naturally expect that all solutions are asymptotically flat. 

i.oK 




Figure 4: Plots of (continuous) and 1 + A(3 (dotted) for different values of ^ < 0, with M = f . From 
left to right, A = ~1 (brown), A = —6 (blue), the limiting case A = — 54e~^ (green), and A = —10 
(red). Notice that the green curve has not precisely a single horizon, but rather two: This shows the 
same numerical limitation exposed in Fig. 2. The dashed black line is the standard general relativity 
result, showing the matching asymptotic behaviour. The zeroes of the dashed lines arc the locations of 
the non-metric singularities of type (1), A/3 = 2. There, diverges, is finite, and thus the metric 
components goo and gn both vanish. 



The plots also show the near horizon structure, the relative positions of the horizon, and 
the non-metric singularity, as zeroes of and 1 + Ap (see also Fig. 2). Notice the presence of 
two zeroes of For A small enough, the non-metric singularity lies inside the event horizon 
given by the outmost zero of f^, i.e., at the outmost zero of 1 -|- Ap > 0. We find that the 
event horizon is at a smaller radius than the Schwarzschild one. As we increase A, the radius of 
the event horizon shrinks, whereas that of the non- metric singularity grows. The critical value 
of A at which the singularity becomes naked can be computed by requiring /^(r^M) < 0, and 
one gets |^| < (|)^. Therefore we have a situation similar to the positive A branch with 
the minimal black hole, having this time the mass Mq = (|)^ /2y^'e\A\. For larger values of |^| 
the non- metric singularity Ap = — 1 is exposed, and the two zeroes of keep moving towards 
each other. The two zeroes merge in a similar manner that two Killing horizons of the Reissner- 
Nordstrom metric merge and form a single degenerate horizon when the extremal limit \Q\ ^ M 
is taken. However, for the present case, both the zeros of occur in the region where 1+Ap < 0: 
both the "horizons" defined as the zeroes of are located inside the non-metric singularity. 
The value at which the two zeros of merge and disappear can be computed by requiring the 
minimum of to be non-negative. We find A < — 54M^/e^ ~ — 7.3fM^ (green line). In [3] this 
relation was interpreted to give, for fixed A, the minimal value of the mass for the black hole 
to exist. However we have seen above that well before (at |^| = (|) ~ — 5.7M^), the black 
hole has already ceased to exist because the non- metric singularity at = — 1 is naked. After 
the zeroes of disappear even the true singularity at r = is exposed. 

Summarising with the help of Fig. 2, we see that in both branches we have at fixed M a 
maximal value of |yl| above which the black hole ceases to exist. This means that for a given A 
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there is always a minimal mass that a black hole can have, Mq = \J Aj (2e^) for positive A and 
Mq = (f )^''^ \J for negative A. In the range of A for which the horizon exists, its location 
at fixed M is an increasing function of A. It is smaller than the Schwarzschild value Th = 2M 
for A < ^ and larger for ^ > 0. The spacetime structure in the case of a naked non- metric 
singularity with no horizons present is shown in the conformal diagram in Figure 5. 




Figure 5: Conformal diagram for the vacuum (Q = = A) solution with a naked non-metric singularity 
(green line), originally shown in [3]. The double line, X, denotes null infinity. The flow of time is vertical, 
different from the original diagram given in [3], and the two green lines should be identified. It is seen 
that the naked non-metric singularity (vertical green line) corresponds to a "Cauchy surface" (horizontal 
green line) for the interior region, and accordingly the central singularity (thick dashed line) is achronal 
(non-timelike) . 
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9 Charged black holes 



We now consider the full system of equations (96) in the presence of an electric charge. The 
numerical studies that we present below show that the metric components approach those cor- 
responding to the Reissner-Nordstr0m metric at large distances, and that (5 goes to zero, thus 

> and the asymptotic signature is H . This confirms the asymptotic flatness of the 

perturbative solution presented earlier (see eq. (98)). Thanks to this asymptotic behaviour, 
which conforms to the conditions (78), (79), (80), the parameters M and Q can be interpreted 
as the ADM mass and the electric charge in the standard sense discussed in Section 6.1. 

Decreasing r from the asymptotic regime, the curvature increases and we start seeing depar- 
tures from the Reissner-Nordstr0m solution in general relativity. In particular, 

• The number and the location of the horizons are changed; 

• There appear again non-metric singularities at A(3 = 2 (of type (1)), or A(3 = —1 (of 
type (2)). 

Concerning the non-metric singularities, they now occur at locations r^M 

(A, M, Q) depending 

also on Q. At these points the metric given by (76) is degenerate or diverging, but as seen in 
the uncharged case, this does not mean that the singularities are harmful to the theory itself. 
In this perspective, we point out that the scalar S constructed in (103) is unchanged by the 
presence of an electromagnetic field, thanks to the property A = 0. This is analogue to 
what happens in general relativity, where the presence of an electric charge does not affect the 
Ricci scalar. Therefore S is again finite at the non-metric singularities. In view of this, although 
we have not fully studied the behaviour of other curvature components at the singularities, we 
anticipate that a similar change of coordinates as the one made in [3] would show finiteness of 
the fundamental fields of the theory. 

As in the uncharged case, the solutions in radial coordinate r interrupt at Ap = 2, where 
/?' diverges (see eq. (96b)), whereas they continue after A(5 = — 1. So again the range of our 
plots is r € (0, oo) if no type (1) singularity occurs, whereas it is {r^yi{Al3 = 2),oo) if it does. 
However, unlike the uncharged case, both types of non-metric singularities can occur regardless 
of the sign of A, since now the sign of /3 is not fixed (see e.g., the perturbative solution (98)). 

To get a first feeling of departures from the general relativity case, in Figure 6, we plot 
the numerical solutions of goo for A = M = 1 and various values of the charge, and compare 
them with the Reissner-Nordstr0m solution in general relativity. These plots show the matching 
asymptotic behaviour anticipated above. The behaviour of goo near the horizon, on the other 
hand, becomes rather different in the modified theory. Note that zeroes of goo are not necessarily 
horizons, as they can also correspond to a non- metric singularity of type (2). 

As we vary A, the departures from general relativity change significantly, and we have to 
face rich fauna of solutions. To give a useful overview of the situation, we plot in Figure 7 the 

^'*The analytic perturbative solution (98) allows us to identify M and Q as the ADM total mass and charge. 
However for practical purposes in the numerical code one has to assign the boundary conditions at a finite value 
of the radius, say tq. The larger is ro, the closer the boundary values M and Q are to the ADM charges; but also 
the less precise are the numerical solutions near the horizon. Conversely, a small ro increases the precision of the 
solution near the horizon, but limitates the interpretation of M and Q. After investigating the stability of our 
code (by comparing it with the analytical cases known at Q = 0, see left panel of Fig. 4) we took ro = 50. 



30 




Figure 6: The numerical solutions for M = 1 and A = 1 oi the modified theory (solid lines) plotted 
against the Reissner-Nordstr0m solution of general relativity (dashed lines) for Q = 0.5 (blue), Q = 0.81 
(green), Q — 1 (orange) and Q = 1.5 (red). Note the matching asymptotic behaviour, goo diverges at 
minus infinity at the radius where the type (1) A(3 = 2 singularity occurs (blue line), whereas diverges 
at plus infinity at r = 0, as the Reissner-Nordstr0m solution, when this type (1) singularity is avoided. 
In the latter case, extra zeroes may appear in goo due to crossings of the type (2) Af3 = — 1 singularity. 
The value Q = 0.81 estimates the turning point from one regime to the other. 



locations of horizons and non-metric singularities as functions of Q, for a few significant values 
of A. We also illustrate the global spacetime structures for a few specific cases in Figures 9, 10, 
and 11. 

The initial parameter space in which the numerical integrations are performed is four- 
dimensional: {^o,A,M,Q). However <l>o and M affect the solutions in a trivial way: a non- 
vanishing cosmological constant <I>o can be straightforwardly added using (96c), and when chang- 
ing the boundary condition M, we simply reproduce the same features at appropriately rescaled 
A and Q. Therefore, we can without loss of generality restrict the parameter space to the 
two-dimensional {A,Q). 

Given the plethora of possibilities, it is useful again to discuss separately the two branches 
with A positive or negative. 

9.1 Positive branch 

Looking at Fig. 6 and the top-left panel of Fig. 7, we can divide the branch with A > in two 
regions: 

• Region (i), where the non-metric singularity of type (1) is present, and all plots interrupt 
at its location rNM(^/3 = 2); 

• Region (ii), where the type (1) singularity is absent, and all plots continue to r = 0. 

For given A, the demarcation between the two regions is given by a certain critical value of 
the charge, which we call Qq{A): when below this value we are in region (i), and when above 
we are in region (ii). In region (i), the location r^M of the non-metric singularity is a decreasing 
function of Q. This decreasing behaviour continues until the critical charge, e.g., Qo ^ 0.8 for 
A = 1, where a minimum value of r^M is reached and we pass in region (ii). In region (ii), the 
type (1) singularity is replaced by the type (2), and its location is now an increasing function of 



31 




Q. Numerical analysis shows that the value Qq{A) which marks the turning point from region 
(i) to region (ii) is a slowly growing function of A. For instance Qo(O-l) ~ 0.55, Qo(l) ~ 0-8 
(see the top-left panel of Fig. 7) and Qq{12) ~ 1.2 (the top-right panel). 

Next, let us move our attention to the horizons, drawn in blue. These are zeroes of 
which are simultaneously infinities of gl- At small Q, the location of the event horizon is a 
decreasing function of Q similar to the one of general relativity, r_|_ = M + \J — Q^, which 
is also shown in the plot for comparison. The difference is that at each Q the location of the 
event horizon is at a larger radius than in the case of general relativity. This generalises to 
the charged case what already seen for Q = in the previous section. A characteristic feature 
of the Reissner-Nordstr0m solution is the presence of an inner horizon. The top-left panel of 
Figure 7 shows that in the modified theory there is only a small range of parameters in which 
the inner horizon is present, and this happens within region (ii). Furthermore, notice that the 
value of the charge at which the two horizons merge is slightly larger than the Q = M value of 
the Reissner-Nordstr0m solution. This result suggests the possibility of screening effects of the 
electric charge in this region of parameter space. 

A stronger departure from general relativity is the presence of a third horizon: A close look 
at the top-left-panel of Figure 7 tells that for Q in a narrow range around Q = 1, we have 
three horizons, the event-, inner-, and the third-horizon (see also right panel of Fig. 6 and Fig. 
11 for the global structure of the spacetime). Then, when Q increases, the first two (event 
and inner) horizons merge and cease to exist (as in the Reissner-Nordst0m case), and only the 
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third one continues to exist for larger value of Q — passed the critical value at which the third 
horizon first appears, and is always behind the type (2) non-metric singularity. The behaviour 
of horizons around Q = 1 and the appearance of the third horizon can be much clearly seen 
for smaller values of A than A = 1. In Figure 8, we show plots similar to Figure 7 but with 
A = 0.1. The locations of both the horizon and non-metric singularity grow with Q, but the 
singularity grows faster. This suggests that there is a region of parameter space in the theory 
where we can arbitrary overcharge a black hole without ever exposing the r = singularity, i.e., 
the central singularities are always behind non-metric singularities (or otherwise covered by an 
event horizon). 

When we increase A, the value Qq{A) at which we pass from region (i) to region (ii) slowly 
grows. Also, the location itself of the non-metric singularities grows with A. Specifically, this 
grows faster than the location of the event horizon, a phenomenon already seen in the uncharged 
case that led to the existence of a notion of "minimal" black hole. The same thing happens in 
the charged case and, because the presence of a charge shrinks the radius of the event horizon, 
we can expect that at fixed M the non-metric singularity of type (1) is exposed at smaller values 
of A as we increase the charge. This situation is depicted in the top-right panel of Fig. 7, where 
we consider the value A = 12. This value is smaller than the A ~ 14.78 required - for M = 1 
- to expose the non-metric singularity in the uncharged case (cf. Fig. 2), but can make the 
type (1) singularity naked if enough charge is present (e.g., Q ~ 0.8 for A = 12, see Figure 7). 

For such a large A (still the top-right panel), we also notice the complete disappearance of 
the event horizon, swallowed up by the type (1) non-metric singularity. This does not exclude 
the possibility than an horizon is still present beyond the singularity, but we do not have access 
to that region. If we continue to increase A, the event horizon disappears for smaller and smaller 
values of Q, until the value A = 2e^M^ found in the previous Section. This is the value for 
which the event horizon disappears also in the uncharged case. Note however that a horizon is 
still present in region (ii), behind the non-metric singularity of type (2). Therefore even at such 
large A, we can still have a horizon that covers the true singularity at the centre, if the charge 
is high enough. 

r 

2.5r 




Figure 8: Locations of the event horizons, inner horizon, and the appearance of third horizon at A = 0.1 
and M = 1. 

In Figures 9, 10 and 11 we report the Penrose diagrams of different configurations for which 
the solution is complete, i.e. r S (0,oo). 
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Figure 9: Conformal diagram for the solution with A = 1 and M ~ 1, Q > 1.1 or A ~ 12, M — 1, 
Q > 1.2. See the top two panels of Fig. 7. The double line, T, denotes null infinity. There is a naked non- 
metric singularity of type (2) (green line), and inside it, a horizon (blue line) that covers a true singularity 
at the centre (thick dashed line). Here, the flow of time is taken to be vertical. As a consequence, in this 
figure, the two green lines — one vertical and the other horizontal — should be identified, as they denote 
the same non-metric singularity. Covered by a horizon (blue line) inside the non-metric singularity, the 
central singularity at r = is non-spacelike. The full extension of the spacetime is obtained by gluing 
together copies of the fundamental block depicted here, at blue solid line. 




Figure 10: Conformal diagram for the charged black hole solution with A ~ 1 and AI = 1, Q = 0.9. 
See the top-left panel of Fig. 7. The double line denotes null infinity. There is a type (2) non-metric 
singularity (green line) hidden by the event horizon (blue line). Inside the non-metric singularity is the 
central singularity (thick dashed line). The flow of time is vertical, and the two green lines — one vertical 
and the other horizontal — are identified as indicated. The central singularity inside the non-metric 
singularity is thus non-spacelike. The full extension of the spacetime is obtained by gluing together 
copies of the fundamental block depicted here, at blue solid line. 
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Figure 11: Conformal diagram for the charged black hole solution with A = 0.1 and M = 1, Q = 0.9. 
See Fig. 8. The double line denotes null infinity. There is a type (2) non-metric singularity (green 
line) inside the inner horizon, and the third-horizon is located inside the non-metric singularity. The 
flow of time is vertical, and the two green lines — one vertical and the other horizontal — denoting the 
same non-metric singularity are identified. The central singularity (thick dashed line) inside the third- 
horizon is non-spacelike. The full extension of the spacetime is obtained by gluing together copies of the 
fundamental block depicted here, at blue solid line. 
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9.2 Negative branch 

In the negative branch the standard two horizons of the Reissner-Nordstr0m solution are more 
easily distinguishable (see the bottom panels of Figure 7). The effect of a negative A to shrink 
the location of the event horizon seen in the uncharged case generalises to the charged case in 
a way that qualitatively maintains the standard dependence on Q. At the same time, the inner 
horizon is pushed to a larger radius. As a consequence, the charge at which the two horizons 
merge and disappear is smaller than the Q = M value of the general relativity case. This hints 
at possible effects of antiscreening of the electric charge in this theory. 

Concerning the appearance of non-metric singularities, we see from the bottom panels of 
Fig. 7 that 

• The solutions always possess the type (1) singularity; 

• The type (2) singularities can also occur in the solutions, for Q smaller than a certain 
value, Qi{A). 

Comparing the two bottom panels of Fig. 7, we see that the value Qi{A) at which the two 
singularities of type (2) merge and disappear grows with |^|. The simultaneous presence of the 
two singularities makes the near horizon behaviour differ significantly from the case of general 
relativity: See the blue plots in Fig. 12 which have A = —1, M = 1 and Q = 0.5. 




Figure 12: The numerical solutions for M — 1 and A ~ —1 oi the modified theory (solid lines) plotted 
against the Reissner-Nordstr0m solution of general relativity (dashed lines) for Q = 0.5 (blue), Q = 1 
(green) and Q = 1.5 (red). Notice the matching asymptotic behaviour, goo always diverges at the 
radius where the A(3 = 2 singularity occurs. For (3 = 1 (green line) there is no horizon, unlike in the 
Reissner-Nordstr0m solution. For Q = 0.5 (blue fines) the singularity at Af3 = —1 is crossed twice, and 
consequently the plot for 1700 shows two additional zeros before the divergence, which is now at minus 
infinity. 

Let us look at the bottom-right panel of Fig. 7. We see that when increasing |^|, the 
two horizons approach each other, and the extremal case happens at a smaller and smaller 
charge. The location of the outmost horizon still maintains the qualitative dependence on Q 
of the Reissner-Nordstr0m solution, but the location of the inner one has now a very deformed 
dependence. The two horizons disappear completely at the value A = — 54M^/e^ computed in 
the uncharged case, and as already seen then, the type (2) singularity becomes exposed before 
the two horizons disappear: the bottom-right panel shows that at ^ = — 7 the horizons are still 
present but inside the singularity of type (2). 
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10 Summary 



We have studied spherically symmetric black holes in the framework of Krasnov's modified self- 
dual gravity [2], in which a Weyl curvature dependent functional is added to the Plebansky 
action of general relativity. The remarkable novelty of this modification is that there appear no 
extra physical degrees of freedom for gravity. This is a significant difference from other modified 
theories of gravity proposed in the literature, and it requires the use of a self-dual 2-form B as the 
fundamental field: The metric describing the geometry of spacetime is only a derived quantity. 
The derived metric can be obtained using the self-duality of B, and the field equations it satisfies 
are still second order as in general relativity. At least for the case of spherical symmetry, the 
field equations can be casted in a way which is significantly similar to general relativity, with the 
effect of the modification clearly distinguishable. The modification increases the non-linearity 
of the equations, due to the extra feedback effect coming from the additional Weyl curvature 
dependent — otherwise "cosmological constant" — term. 

The first result of the paper is to show how to couple electromagnetism to the theory. This 
uses the coupling to Plebahski's theory developed in [10], which is shown to apply straightfor- 
wardly also to the modified theory, as anticipated by Krasnov [2]. We have then proved the 
staticity of the electro-vacuum spherically symmetric solutions for any form of $ (Birkhoff 's 
theorem, Section 5.6). Under the assumption of spherical symmetry and accordingly with the 
staticity result, we have considered spherically symmetric, static black hole spacetimes for the 
two simplest profiles of i.e. the linear modification <1> = —A/3 + a/5(r) and the quadratic 
modification $ = —A/3 + Af3'^{r)/2, with a and A parameters and /3(r) the 'magnitude' of Weyl 
curvature of our spherically symmetric spacetimes. In both cases, the structures near the hori- 
zon and near the central singularity significantly differs from general relativity, and the notion 
of extremal black hole changes with the parameters. We summarise our main results below. 

In the linearly modified case, the solutions do not approach asymptotically flat geometries 
at large distances even when the cosmological constant $o = —A/3 is set to zero, as we have 
briefly discussed in Section 6.1, and hence the solutions have no direct physical relevance. Since 
the linear modiflcation fails to satisfy a certain notion of analyticity, the question arises if there 
is a relation between the latter and the possibility of asymptotic flatness. Such a relation 
would be useful to restrict the class of <^'s with physically relevant solutions. The linear case is 
nonetheless worth examining since it can be easily solved analytically. The solutions obtained 
in Section 6.2 explicitly show how the Q-dependence of the locations of the horizons is affected 
by the parameter a. 

In the quadratically modified case, the non-linearity of the equations of motion is more in- 
volved, and we have not been able to find analytic solutions for Q ^ 0. We resorted to numerical 
computations to study the charged black holes. For this purpose, we first considered a pertur- 
bative approach to the field equations. The approximate solutions found in this way are valid at 
large distances, and allowed us to show that the solutions have the standard Reissner-Nordstr0m 
asymptotic behaviour, and thus the parameters (M, Q) appeared there can be interpreted in the 
standard sense, the ADM mass and the total charge, respectively, in accordance with the argu- 
ments of Section 6.1. 

The numerical methods allowed us also to further investigate the quadratic modification case 
with Q = = A, whose vacuum black hole solutions had already been analytically found in [3]. 
We studied the location of the event horizon, and showed that a black hole of mass M can only 
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exists for A G ( — (|)'^e~^M^, 2e^M^), where the lower bound differs from the one proposed in 
[3]. For values of A outside this interval, a non-metric singularity is naked outside the horizon; 
it is of type (1) if ^ > 0, and of type (2) if ^ < 0. Concerning the nature of the singularity, 
it was shown in [3] that a change of coordinates can be found in which the fundamental fields 
are finite. Therefore this is not a physical singularity for the theory, and the fundamental fields 
can be continued beyond it. Here we have introduced also a scalar quantity, analogue of the 
Ricci scalar, and showed explicitly that it is finite at the non- metric singularities. The vacuum 
Q = = A case results are summarised as follows. 

1. The location of the horizon, rii{A, M), at fixed M is pushed at larger /smaller values than 
in GR according to the sign of A. 

2. The internal structure becomes more intricate: the true singularity is surrounded by the 
non-metric singularity and an internal horizon. The non- metric singularity is of type (1), 
Ap = 2, for yl > and type (2), A/3 = —1, for yl < 0. The new internal horizon can be 
seen explicitly for A < 0, but for A > one has to change branch of the solution after the 
singularity to see it, or go to /3 as radial coordinate, as it is done in [3]. 

3. For too large, the non-metric singularity is exposed outside the event horizon. Therefore 
in a given modified action with fixed A, there is a minimal mass that a black hole can 
have. Beyond this limit, distant observers would not see the event horizon, but a non- 
metric singularity first. 

We next summarise the charged case. Our main focus has been on which type of (e.g., (1) 
Af3 = 2 or (2)A/3 = —1) non-metric singularities appear and whether those singularities are 
covered by the event horizon. Given a fixed M, the answer depends upon the values of the 
parameters A and Q. 

1. As in GR, one can find an event horizon as well as an inner horizon. The location 
rii{A,Q, M) of the event horizon for fixed M and Q is pushed at larger or smaller val- 
ues than in GR according to the sign of A. For the inner horizon, we have the opposite 
behaviour. As a consequence, the extremal CQjSG is how 3jt 3i \3iT ger or smaller value than 
Q = M as in GR. This suggests that in the theory, electromagnetic charges might experi- 
ence the screening or antiscreening effects according to the sign of A, as we have discussed 
in Section 9. The physical implication of the screening effects, however, remains an open 
issue. 

2. The true singularity at r = is surrounded by non-metric singularities, and (if exists) 
inner horizons. The outmost non-metric singularity is: 

• type (1) at small Q, then (2) at larger Q for ^4 > 

• typs (2) at small Q, then (1) at larger Q for ^4 < 

3. Increasing \A\ too much, the non-metric singularity is exposed outside the event horizon. 
This happens earlier and earlier, as we pour a larger charge into the hole. Therefore in a 
given modified action with fixed A, there is a minimal mass that a black hole of charge Q 
can have, and this value increases as we increase Q. 
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4. A region of parameter space exists, with ^ > 0, where the type (2) singularity is naked, but 
inside it there always exists a horizon which hides the singularity at r = for arbitrarily 
large values of the charge. 

Note that there could be further structure behind all singularities of type (1), which we do 
not attempt to access numerically. This would require a change of coordinates, or the knowledge 
of the analytic solution as it was done for the vacuum solution in [3]. (Compare with Section 
8 for the positive branch of the uncharged case.) It would certainly be interesting to do so, 
changing the radial coordinate as in the vacuum case. However, given the limitations of our 
understanding of analytic solutions, we prefer to restrain from doing so, and focus instead the 
attention of the reader on the near horizon structure. 

Our results also bring some insight into how it is possible to restrict the class of physically 
relevant types of <I>. We showed that spherically symmetric black holes (charged or not) are not 
asymptotically flat in some type of <1> non-analytic at the origin, and one is tempted to conjecture 
that this generalises to any non-analytic choice. Although charged black holes are not relevant 
in astrophysical context, the results obtained show that a deeper understanding of the modified 
theory, as well as further useful restrictions on $ can be found by looking at matter coupling. 

Finally, we would like to recall that given M, there are parameter regions of {A, Q) in which 
both types of non-metric singularities can be naked. In this case, distant observers would see the 
breakdown of the standard spacetime picture based on the metric description outside the horizon. 
Nevertheless, from the view point of this modified self-dual gravity, there is nothing to be afraid 
of, since the non-metric singularities are harmless for the fundamental fields, even when naked. 
Therefore, in this modified self-dual gravity, the notion of the cosmic censorship itself may also 
need be modified. In this regard, we should also recall that a non-metric singularity is specified 
as a (r = const.) hypersurface, across which the time and spatial (radial)-coordinates change 
their spacetime role. This, in particular, implies that as illustrated by the conformal diagram 
in Figure 5 for the simplest vacuum case, a true singularity at the centre (r = 0) covered by a 
naked non-metric singularity is achronal (non-timelike), unless covered by any other horizons. 
We also note from Figure 5 that the hypersurface of a naked non-metric singularity plays the 
role of the (partial) Cauchy surface for the interior region, and therefore once a causal curve goes 
across a non-metric singularity and enters the interior region, then it has to fall into the central 
singularity and never comes back to the region outside the non-metric singularity. Therefore, in 
this sense the hypersurface of a non-metric singularity may be viewed as a one-way membrane, 
or a type of "horizon." This view point may appear to be in favour of the cosmic censorship 
in the sense that the true singularities of the theory appear to always be covered by either the 
event horizon, or a non-metric singular horizon. Note that in the present context, one cannot 
tell which horizon, future or past, a non-metric singularity corresponds to, due to the time- 
symmetric nature of static spacetimes. However, the situation would be clear when dynamical 
(non-stationary) cases are considered. It would be interesting to consider a sensible definition 
of a Cauchy problem for the evolution in such a spacetime and the notion of cosmic censorship, 
within the framework of the modified self-dual gravity theory. 
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